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Abstract. This paper investigates several global rigidity issues for polyhedral surfaces 
including inversive distance circle packings. Inversive distance circle packings are poly- 
hedral surfaces introduced by P. Bowers and K. Stephenson in [2] as a generalization of 
Andreev-Thurston's circle packing. They conjectured that inversive distance circle pack- 
ings are rigid. Using a recent work of R. Guo [9] on variational principle associated to 
the inversive distance circle packing, we prove rigidity conjecture of Bowers-Stephenson 
in this paper. We also show that each polyhedral metric on a triangulated surface is de- 
termined by various discrete curvatures introduced in [12., verifying a conjecture in |12j . 
As a consequence, we show that the discrete Laplacian operator determines a Euclidean 
polyhedral metric up to scaling. 



1. Introduction 

1.1. This is a continuation of the study of polyhedral surfaces |12j . |13j . The paper 
focuses on inversive distance circle packings introduced by Bowers and Stephenson and 
several other rigidity issues. Using a recent work of Ren Guo jU], we prove a conjecture of 
Bowers-Stephenson that inversive distance circle packings are rigid. Namely, a Euclidean 
inversive distance circle packing on a compact surface is determined up to scaling by its 
discrete curvature. This generalizes an earlier result of Andreev [TJ and Thurston |17] on 
the rigidity of circle packing with acute intersection angles. In j!2j . using 2-dimensional 
Schlaefli formulas, we introduced two families of discrete curvatures for polyhedral surfaces 
and conjectured that each of one these discrete curvatures determines the polyhedral metric 
(up to scaling in the Euclidean case). We verify this conjecture in the paper. One conse- 
quence is that for a Euclidean or spherical polyhedral metric on a surface, the cotangent 
discrete Laplacian operator determines the metric (up to scaling in the case of Euclidean 
metric). The theorems are proved using variational principles and are based on the work 
of [H] and j!2j . The main idea of the paper comes from reading of @], [Z] and [15j. 
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1.2. Recall that a Euclidean (or spherical or hyperbolic) polyhedral surface is a tri- 
angulated surface with a metric, called a polyhedral metric, so that each triangle in the 
triangulation is isometric to a Euclidean (or spherical or hyperbolic) triangle. To be more 
precise, let E 2 , S 2 and H 2 be the Euclidean, the spherical and the hyperbolic 2-dimensional 
geometries. Suppose (S, T) is a closed triangulated surface so that T is the triangulation, 
E and V are the sets of all edges and vertices. A K 2 (K 2 = E 2 , or S 2 , or H 2 ) polyhedral 
metric on (S, T) is a map / : E — > R so that whenever e,, e^, e k are three edges of a triangle 
in T, then 

l(ei) + l(ej) > l(e k ), 

and if K 2 = S 2 , in addition to the inequalities above, one requires 

l( ei ) + l( ej ) + l(e k ) < 2tt. 

Given I : E — > R satisfying the inequalities above, there is a metric on the surface S, called 
a polyhedral metric, so that the restriction of the metric to each triangle is isometric to 
a triangle in K 2 geometry and the length of each edge e in the metric is /(e). We also 
call I : E — y R the edge length function. For instance, the boundary of a generic convex 
polytope in the 3-dimensional space E 3 , or S 3 or H 3 of constant curvature 0, 1, or —1 is a 
polyhedral surface. The discrete curvature k of a polyhedral surface is a function k : V — > R 
so that k(v ) = 2tt — YlT=i ^ i where O^s are the angles at the vertex v. See figure 1. 

Since the discrete curvature is built from inner angles of triangles, we consider inner 
angles of triangles as the basic unit of measurement of curvature. Using inner angles, we 
introduce three families of curvature like quantities in |12j . The relationships between the 
polyhedral metrics and curvatures are the focus of the study in this paper. 

Definition 1.1. Q12J) Let h e R. Given a K 2 polyhedral metric on (S,T) where K 2 
= E 2 , or S 2 or H 2 , the <ph curvature of a polyhedral metric is the function <p h : E — > R 
sending an edge e to: 

pa pa' 

(1.1) Me)= / sm h (t)dt+ / sm h (t)dt 

Jn/2 A/2 

where a, a' are the inner angles facing the edge e. See figure 1. 

The i/jh curvature of the metric I is the function iph : E — > R sending an edge e to 

b+c-a b'+c'-a' 

(1.2) V>fc(e)= / 2 cos h (t)dt+ [ 2 cos h {t)dt 

Jo Jo 

where b, b', c, d are inner angles adjacent to the edge e and a, a' are the angles facing the 
edge e. See figure 1. 

The curvatures 0o an d ipo were first introduced by I. Rivin [Ri] and G. Leibon [Le] 
respectively. If the surface S = S 2 , then these curvatures are essentially the dihedral 
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Figure 1. 

angles of the associated 3-dimensional hyperbolic polyhedra at edges. The curvature 
0_2(e) = — cot (a) — cot(a') is the discrete (cotangent) Laplacian operator on a polyhedral 
surface derived from the finite element approximation of the smooth Beltrami Laplacian 
on Riemannian manifolds. 

One of the remarkable theorems proved by Rivin |15j is that a Euclidean polyhedral 
metric on a triangulated surface is determined up to scaling by its <po discrete curvature. 
In particular, he proved that an ideal convex hyperbolic polyhedron is determined up to 
isometry by its dihedral angles. 

We prove, 

Theorem 1.2. Let (S,T) be a closed triangulated connected surface. Then for any h G R, 

(1) a Euclidean polyhedral metric on (S, T) is determined up to isometry and scaling by 
its 4>h curvature. 

(2) a spherical polyhedral metric on (S,T) is determined up to isometry by its <ph cur- 
vature. 

(3) a hyperbolic polyhedral surface is determined up to isometry by its iph curvature. 

We remark that theorem 1.2(1) for h = was aforementioned Rivin's theorem. However, 
our proof of Rivin's theorem is different from that in |15j and we use the variational principle 
established by Cohen-Kenyon-Propp [5]. Theorem 1.2(3) for h = was first proved by 
Leibon [11] . Theorem 1.2(2) for h = was proved in [14] and theorem 1.2(2) and (3) for 
h < —1 or h > was proved in [12] . 

Take h = —2 in theorem 1.2, we obtain, 

Corollary 1.3. (1) A connected Euclidean polyhedral surface is determined up to scaling 
by its discrete Laplacian operator. 

(2) A spherical polyhedral surface is determined by its discrete Laplacian operator. 

Note that for a Euclidean polyhedral surface, <ph = iph- There remain two questions on 
whether 0^ curvature determines a hyperbolic polyhedral surface or whether ip^ curvature 
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determines a spherical polyhedral surface. It seems the results may still be true in these 
cases. 



1.3. Inversive distance circle packings are polyhedral metrics on a triangulated surface 
introduced by Bowers and Stephenson in [2]. An expansion of the discussion of [2] is in 
[3J. See also |16j . They are generalizations of Andreev and Thurston's circle packings. 
Unlike the case of Andreev and Thurston where adjacent circles are intersecting, Bowers 
and Stephenson allow adjacent circles to be disjoint and measure their relative positions 
by the inversive distance. As observed in [2], this relaxation of intersection condition is 
very useful for practical applications of circle packing to many fields, including medical 
imaging and computer graphics. Based on extensive numerical evidences, they conjectured 
the rigidity and convergence of inversive distance circle packings in [2]. Our result shows 
that Bowers-Stephenson's rigidity conjecture holds. The proof is based on a recent work of 
Ren Guo [S] which established a variational principle for inversive distance circle packings. 
A very nice geometric interpretation of the variational principle was given in [5]. 

We begin with a brief recall of the inversive distance in Euclidean, hyperbolic and 
spherical geometries. See [3] for a more detailed discussion. Let K 2 be E 2 , or H 2 or S 2 . 
Given two circles C\,C?, in K 2 centered at t>i,t>2 of radii r\ and r 2 so that t>i,t> 2 are of 
distance / apart, the inversive distance / = I(C\,C2) between the circles is given by 

(1.3) I= P-r 2 -r 2 



2rir 2 

in the Euclidean plane, 

,j ^. j. cosh(Z) — cosh(ri) cosh(r 2 ) 

sinh(ri) sinh(r2) 

in the hyperbolic plane and 



cos(/) - cos(ri) cos(r 2 ) 
sm(ri) sm(r 2 ) 

in the 2-sphere. See [9] for more details on (1.4) and (1.5). If one considers E 2 , H 2 and S 2 
as appeared in the infinity of the hyperbolic 3-space H 3 , then G\ and C 2 are the boundary 
of two totally geodesic hyperp lanes D\ and .D 2 . The inversive distance I is essentially the 
hyperbolic distance (or the intersection angle) between D\ and _D 2 . In particular, for the 
Euclidean plane E 2 , the inversive distance J(Ci,C 2 ) is invariant under the inversion and 
hence the name. 

Bowers and Stephenson's construction of an inversive distance circle packing with pre- 
scribed inversive distance on a triangulated surface (S, T) is as follows. Fix once and for 
all a vector I G [ — 1, oo) E , called the inversive distance. 
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In the Euclidean case, for any r G R> 0) called the radius vector, define the edge length 
function / G Rf by the formula 

(1.6) /(e) = y/r~{v) 2 + r(u) 2 + 2r(v)r(u)I(e) 

where the end points of the edge e is {u, v}. If /(e)'s satisfy the triangular inequalities that 

(1.7) l(e i ) + l(e j )>l(e k ) 

for three edges e*, ej, e k of each triangle in T, then the length function I : E — >■ R sending e 
to /(e) defines a Euclidean polyhedral metric on (S,T) called the inversive distance circle 
packing with inversive distance /(e) at edge e. Note that if 1(e) G [0, 1] for all e, the 
polyhedral metric is the circle packing investigated by Andreev and Thurston where the 
intersection angle between two circles at the end points of an edge is arccos(/(e)). 

In the hyperbolic geometry, one uses 

(1.8) /(e) = cosh _1 (cosh(r(i>)) cosh(r(w)) + /(e) sinh(r(t> )) sinh(r(w)) 

as the length of an edge. If (1.7) holds, then the lengths 1(e) 's define a hyperbolic inversive 
distance circle packing with inversive distances / on (S, T). The spherical inversive distance 
circle packing is defined similarly with additional condition on /(e) 's that 

1(d) + l( ej ) + l(e k ) < 2tt 

for each triangle with edges e iy ej, e k . 

The geometric meaning of these polyhedral metrics is the following. In each metric, if 
one draws a circle of radius r(v) at each vertex v, then inversive distance of two circles at 
the end points of an edge e is the given number /(e). 

Our result which solves Bowers-Stephenson's rigidity conjecture is the following. 

Theorem 1.4. Given a closed triangulated connected surface (S,T) with the set of edges 
E and I G R> considered as the inversive distance, 

(1) a hyperbolic inversive distance circle packing metric on (S,T) of inversive distance 
I is determined by its discrete curvature k : V — > R. 

(2) an Euclidean inversive distance circle packing metric on (S,T) of inversive distance 
I is determined by its discrete curvature k : V — >■ R up to scaling. 

Note that for / G [0, 1] E , the above result was Andreev-Thurston's rigidity for circle 
packing with intersection angles between [0,7r/2]. It seems the similar result may be true 
for / G [—1, oo) E . 

1.4. The paper is organized as follows. In §2, we prove an extension lemma for angles 
of triangles. We also establish a criterion for extending a locally convex function to convex 
function. In §3, we prove theorem 1.4. Theorem 1.2 is proved in §4. 

The following notations and conventions will be used in the paper. We use R, R>o, 
R>o, R<o to denote the sets of all real numbers, positive real numbers, non- negative real 
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numbers, and negative real numbers respectively. If X is a set, R x = {/ : X — > R} is the 
vector space of all functions on X. If A is a subspace of a topological space X, then the 
closure of A in X is denoted by A. 

We thank Ren Guo for comments and careful reading of the manuscript. 



2.1. Continuous extension by constants. 

Definition 2.1. Suppose A is a subspace of a topological space X and / : A — >■ Y is 
continuous. If there exists a continuous function F : X — > Y so that F\a — f and F is a 
constant function on each connected component of X — A, then we say / can be extended 
continuously by constant functions to X. 

Note that if each connected component of X — A intersects the closure of A, then the 
extension function F is uniquely determined by /. 

The key observation of the paper is the following simple lemma. 

Lemma 2.2. Suppose A is a triangle in the Euclidean plane E 2 7 or the hyperbolic plane 
H 2 7 or the 2-sphere S 2 so that its edge lengths are l±, l 2 , h and its inner angles are 9±, 9 2 , 9 3 . 
Assume that 9; L 's angle is opposite to the edge of length U for each i. Consider Q { = 9i(l) as 
a function of I = (li,h,h)- 

(1) If A is Euclidean or hyperbolic, the angle function 9i defined on Q = {(li,l 2 ,h) £ 
R 3 |^i + k > I3J2 + I3 > h,h + h > h} can be extended continuously by constant 
functions to a function 9>i on R> - 

(2) If A is spherical, the angle function 9i defined on Q = {(h,l2,h) £ R 3 |^i + h > 
h,h + h > h,h + h > hji + h + h < 2n} can be extended continuously by constant 
functions to a function 9{ on (0,7r) 3 . 

We call the set Q in the lemma the natural domain of the length vectors. 

Proof. In the case (1), the extension function 9i of 9i is given by 9i — tt when U > 
and 9i = when lj > k + h- It remains to verify the continuity of 9i on R^q- It is based 
on the cosine law. Given a point L = L 2 , L 3 ) in the boundary Cl — fl of Q inside R^q, 
we may assume without loss of generality that L 1 = L 2 + L 3 . The continuity of 9^ follows 
from 



2. Convex Extension of Locally Convex Functions 



n0i(Z)=7r, lim^(/) = 0, j 

L I — vL 



7 = 2,3. 



Indeed, the cosine law says, in the case of A C E 2 , that 



(2.1) 
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One sees easily that when / tends to L, then the right-hand-side of (2.1) tends to 1 if 
i=2,3 and —1 if % = 1. This verifies the continuity in the Euclidean case. In the hyperbolic 
case, the cosine law says 

(2.2) oos^) = C ° Sh(g C " M4) " " SHk) - 

K J K J sinh(^) sinh(Z fc ) 

Thus one sees that as / tends to L = (Li, L 2 , L3) with Lj > 0, the right-hand-side of (2.2) 
tends to 1 if % — 2, 3 and to — 1 if % — 1. Thus is continuous. 

To see (2), recall that the cosine law for spherical triangle says 

cos(Zj) -cos(fj)cos(f fc ) 

K z - 6 ) cos{t>i) = , , . 

sm(y sin(/ fe ) 

If Z tends to L where L x = L 2 + L 3 with Lj G (0, 7r), then linu_>£ cos(#i) = —1 and 
lim^i cos(ftj) = 1 when % — 2, 3. On the other hand, if L x + L 2 + L 3 = 2n for Lj G (0, 7r), 
then the cosine law implies that lim^i cos(^) = —1 for all i, i.e., all inner angles are 7r in 
this case. Thus by setting the extended function 9^ in (0, 7r) 3 to be Q^l) — n if > lj + l^, 
9i(l) = if lj > k + h, and 8^(1) — ir if ^ + lj + l k > 7r, ( {i,j, k] = {1, 2, 3}), we see that 
6i is continuous. 

□ 

2.2. Continuous extension of 1-forms and of locally convex functions. We 

establish some simple facts on extending closed 1-forms and locally convex functions to 
convex functions in this subsection. 

Definition 2.3. A differential 1-form w = Y17=i a i{ x )d- x i m an open set U C R ra is said to 
be continuous if each a^x) is a continuous function on U . A continuous 1-form w is called 
closed if J 9t w = for each Euclidean triangle r in U . 

By the standard approximation theory, if w is closed and 7 is a piecewise smooth 
null homologous loop in U, then J^w = 0. If U is simply connected, then the integral 
F(x) = f^w is well defined, independent of the choice of piecewise smooth paths in U 
from a to x. The function F(x) is C 1 -smooth so that dF(x)/dxi = aj(x). 

Proposition 2.4. Suppose X is an open set in R n and A C X is an open subset bounded 
by a smooth (n-1)- dimensional submanifold in X. If w = J27=i a i( x )^ x i ^ s a continuous 
1-form on X so that w\a and w\ X -A are closed where A is the closure of A in X, then w 
is closed in X. 

Proof. Since closedness is a local property and is invariant under smooth change of 
coordinates in X, we may assume that X = R n and A = {(xi, ...,x n ) G R n |x n > 0}. Take a 
Euclidean triangle r C X. To verify J di _ w = 0, we may assume that r is not in A or X — A 
since otherwise J dr w = follows from the assumption and the standard approximation 
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theorem. In the remaining case, r intersects both A and X — A. The plane x n = cuts the 
triangle r into a triangle 71 and a quadrilateral 72 so that 71 and 7 2 are in the closure of 
A and X — A. We can express, in the singular chain level, dr = cfyi + 972. By definition, 

f(hi w = oior each *■ Thus far w = Ja 71 w + Ja 72 w = °- D 

A real analytic codimension- 1 submanifold Y in an open set X in R n is a smooth 
submanifold so that locally F is defined by k(x) = for a non-constant real analytic 
function k. Note that if L is a (compact) line segment in X, then either L C F or L fl Y is 
a finite set. This is due to the fact that a non-constant real analytic function on an open 
interval has isolated zeros. 

Recall that a function / defined on a convex set X C R™ is called convex if for all 
p,q E X and all £ G [0, 1], tf{p) + (1 — t)f(q) > f(tp + (1 — t)q). It is called strictly convex 
if for all p ^ q in X and all t G (0, 1), tf(p) + (1 - > f(tp + (1 - t)g). A function / 

defined in an open set U C R n is said to be locally convex (or locally strictly convex ) if it 
is convex (or strictly convex) in a convex neighborhood of each point. 

Proposition 2.5. Suppose X C R n is an open convex set and A C X is an open subset of 
X bounded by a codimension- 1 real analytic submanifold in X. If w = Y^=i a i( x )dxi is a 
continuous closed 1-form on X so that F{x) = f*w is locally convex in A and in X — A, 
then F(x) is convex in X . 

Proof. Since X is simply connected, the function F is well defined. To verify convexity, 
take p,q G X and consider f(t) = F(tp + (1 — t)q) for t G [0, 1]. It suffices to show that 
f(t) is convex in t. Since F is C 1 -smooth, / is C 1 -smooth. Let OA = A — A and L be the 
line segment from p to q. Since dA is real analytic, either L intersects dA in a finite set 
of points, or L is in dA. In the first case, let = t < t\ < ...,£„ = 1 be the partition of 
[0, 1] so that the line segment tp + (1 — t)q for t G (t i: t i+1 ) is either in A or in X — A. By 
definition, f(t) is convex in [tj,t i+1 ], i.e., f'(t) is increasing in [tj,tj +1 ] for % = 0, ...,n — 1. 
Since f'(t) is continuous in [0,1], this implies that f'(t) is increasing in [0,1], i.e., f(t) is 
convex in [0, 1]. In the second case that L C dA, we take two sequences of points p m and 
q m converging to p and q respectively in X so that p m , q m are not in dA. Then by the case 
just proved, the functions f m (t) = F(tp m + (1 — t)q m ) are convex in t. Furthermore, f m 
converges to /. Thus / is convex. □ 

Corollary 2.6. Suppose X C R n is an open convex set and A C X is an open subset 
of X bounded by a real analytic codimension- 1 submanifold in X. If w = J2™ =1 ai(x)dxi 
is a continuous closed 1-form on A so that F(x) = f x w is locally convex on A and each 
a,i can be extended continuously to X by constant functions to a function a~i on X , then 
F(x) = Y17=i didxi is a ^-smooth convex function on X extending F. 

We remark that the real analytic assumption in the proposition 2.5 can be relaxed to 
C 2 smooth. 
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3. A Proof of Bowers-Stephenson's Rigidity Conjecture 



We begin by recalling Guo's work on a variational principle associated to inversive dis- 
tance circle packings and then prove theorem 1.4. We will work in Euclidean and hyperbolic 
geometries only. 

3.1. Guo's variational principle for inversive distance circle packing. Suppose 
A is a triangle with vertices Vi,v 2 ,v 3 and edges = ViVj, i ^ j. Fix once and for all an 
inversive distance G [0, oo) at each edge e^. Then for each assignment of positive number 
Ti at Vi for z = l,2, 3, let 



for hyperbolic geometry where {i,j, k} = {1, 2, 3}. 

Let Q = {(xi, x 2 , £3) G R 3 |xj + Xj > Xk, k} = {1, 2, 3}}. If (/1, l 2 , h) is in Q, then 
we construct a Euclidean triangle A with length Ik of given by (3.1) and a hyperbolic 
triangle, still denoted by A, with length 4 of ey given by (3.2). Suppose the angle of the 
triangle at is Q{ and consider 9{ as a function of (?"!, r 2 , r 3 ). Guo proved the following 
theorem in [9]. 

Theorem 3.1. (Guo $\) Fix any (I 12 , I 23 , hi) G [0, oo) 3 . 

(1) For Euclidean triangles, let Ui = lnr i; then the differential 1-form w = Y^=i^ u i 
is closed in the open subset 0/R 3 where it is defined. The integral F(u) = f^w is a locally 
concave function in u = (wi,^,"^) and is strictly locally concave in u\ + u 2 + M3 = 0. 
Furthermore, if c G R and F(u) is defined, then F{u + (c, c, c)) = F{u). 

(2) For hyperbolic triangles, let = ln(tanh(rj/2)), then the differential 1-form w = 
^2i=i @idui is closed in the open subset o/R 3 where it is defined. Furthermore, the integral 
F(u) = J", j X1 \W is a strictly locally concave function in u = (ux, u 2 , M3). 

It is also proved in [9 J that the open sets where the 1-forms w are defined in theorem 3.1 
are connected and simply connected. Theorem 3.1 is a generalization of an earlier result 
obtained in [6] . Guo proved a local and infinitesimal rigidity theorem for inversive distance 
circle packing using theorem 3.1. It says that a Euclidean inversive distance circle packing 
is locally determined, up to scaling, by the discrete curvature of the underlying polyhedral 
surface. He also proved the local and infinitesimal rigidity for hyperbolic inversive distance 
circle packings. 

3.2. Concave extension of Guo's action functional. Our main observation is that 
Guo's differential 1-forms w = Yll=i ®idui can be extended to a closed 1-form on R 3 in the 
Euclidean case and on R^q in the hyperbolic case so that the integrations of the extended 
1-forms are still concave. 




Ik = cosh (cosh(rj) cosh(r\,-) + sinh(rj) sinh(rj)) 
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Proposition 3.2. Let w be the 1-forms defined in theorem 3.1. 

(a) In the case of Euclidean triangles, the 1-form w can be extended to a continuous 
closed 1-form w on R 3 so that the integration F(u) = J " w is a C 1 -smooth concave function. 

(b) In the case of hyperbolic triangles, the 1-form w can be extended to a continuous 
closed 1-form w on R 3 so that the integration F(u) = /"^ is a C 1 -smooth concave 
function. 

We begin by focusing the 1-forms in its radius coordinate r = (r 1 ,r 2 ,r 3 ) G R> - In 
this case, the 1-forms are given by w = J2i=i ^T 1 anc ^ w = Yli=i ®i sinter ) • ^he l-f° rm w 
is defined on the open set U of R 3 where 

(3.3) U = {(n, r 2 , r 3 ) G R 3 >0 \h + lj > h, k} = {1, 2, 3}}, 

where U = li(ri,r 2 ,r 3 ) is defined on R 3 . (Note that for hyperbolic and Euclidean geome- 
tries, the sets U are different due to (3.1) and (3.2)). The extension of the 1-form w is the 
natural one. Namely, we replace 0t in w by d; t appeared in lemma 2.1. Thus the extended 
1-form is w = 8i— or w = ^-^K- 

It remains to show that w is continuous and closed in R^q so that its pull back to the 
■u-coordinate has a concave integration. To this end, we prove, 

Lemma 3.3. Let U be the closure of U in R> . Then, 

(1) 0i is a constant function on each connected component ofU — U, and 

(2) for each connected component V o/R> — U , the intersection VD U is a connected 
component of U — U . 

PROOF. By (3.3), the boundary dU — U - U is given by U 3 =1 <9j£7 where diU = 
{(ri,r 2 ,r 3 ) G Rl \k = h + h,{j,k} = {1,2,3} - {i}}. Furthermore, R 3 - U = U 3 =1 V^ 
where V t = {(n, r 2 , r 3 ) G R 3 ^ > ^ + Z fc , {j, k} = {1, 2, 3} - {i}}. 

First, we note that if < 1, then dkU = and = 0. Indeed, if Iij < 1, then by 
(3.1) and (3.2), 

^ < n + Tj. 

But due to I ab > 0, (3.1) and (3.2), r,- < ^ and r, < /j. Therefore, l k < k + lj. This implies 
that <9 fc £7 = and V k = 0. 

Next $17 n o^t/ = and ^ n ^ = for i ^ j. Indeed, if r G diU H o^t/ or r G n 
then li > lj + Ik and > /j + Ik- Thus = 0. But Ik > > 0. 

We claim that if > 1, then both \4 and c^f/ are non-empty and connected. Assume 
the claim, then the lemma follows. Indeed, since l s > for all indices s, it follows, by 
lemma 2.1, that 0, is either or ir in d s U, i.e., (1) holds. Next, V^'s are the connected 
components of R 3 — U so that V s fl U — d s U. Thus (2) holds. 

To see the claim, it suffices to show that there is a smooth function /(r^r^) defined on 
R 3 so that its graph is d k U and V k = {(ri,r 2 ,r 3 ) G R> |0 < r 3 < f(r 1 ,r 2 )}. 
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To this end, consider the equation 
(3.4) h = h + lj, 

and let the right- hand- side of (3.4) be g(r k ,ri,rj). We will deal with the Euclidean and 
hyperbolic geometry separately. 

CASE 1 Euclidean triangles. In this case, the function g(r k ,ri,rj) is given by 



(3.5) g(r k , r h r,) = ^Jr 2 k + r 2 + 2I kj r k rj + yj rf + r\ + 2/ iA .r;r fc 

Evidently, for a fixed (r i; rj) G R>o 5 g{ r k, r i, r j) is a strictly increasing function of r k G R>o 



so that g(0, r h rj) = r^+rj < ^Jr 2 + r 2 + 21^^ j (due to Uj > 1) and lim^^ g(r k , n, r k ) = 
oo. By the mean-value theorem, there exists a unique positive number /(r^r^) so that 



g(f(ri,rj),ri,rj) = yrf + r 2 + 2^^/^ = The smoothness of f(r, h rj) follows from the 
implicit function theorem applied to (3.4). Indeed, 

dg = r k + 2I k fj r k + 2I ik rj ^ Q 

<9r fe /i J, 

Thus, f(ri,rj) is smooth. 

This shows is the graph of the smooth function / defined on R> > i- e -> 

d k U = {(ri,r 2 ,r 3 ) G R| |r fc = /(r^r,)}. 

Thus it is connected. Since g(r k ,ri,rj) is an increasing function of r k , V k — {r G -R>o|0 < 
< /(?"*, ?"_?)> {m} = {1,2,3} - {k}}. Thus \4 is connected. 

CASE 2 hyperbolic triangles. By the same argument as in case 1, it suffices to show 
the same properties established in case 1 hold for g(r k ,ri,rj) given by 
(3.6) 

cosh 1 (cosh(r i ) cosh(r k )+I ik sinh(rj) sinh(r fc ))+cosh 1 (cosh(r fc ) cosh(rj)+I k j sinh(r fc ) sinh(rj)). 

Fix (ri,rj) G R> - Then the function g(r k ,ri,rj) is clearly strictly increasing in r k G R>o 
so that linir^oo g(r k , r^ rj) = oo and due to > 1, 

g(0,ri,rj) = r, + rj 
= cosh _1 (cosh(rj + r,)) 
= cosh _1 (cosh(rj) cosh(rj) + sinh(rj) sinh(r,-)) 
< cosh _1 (cosh(rj) cosh(rj) + 1^ sinh(rj) sinh(rj)) = l k . 

By the mean value theorem, there exists a unique positive number /(r^r^) so that 
g(f(ri,rj),ri,rj) = l k . The smoothness of f(ri,rj) follows form the implicit function theo- 
rem that 

dg cosh(rj) sinh(rfc) + I ik sinh(rj) cosh(rfc) 

d r k a/ (cosh(rj) cosh(r fe ) + I ik sinh(rj) sinh(r fc )) 2 - 1 
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cosh(rj) sinh(r fe ) + Ij k sinh(rj) cosh(r fc ) 
a/ (cosh(rj) cosh(r fc ) + Ij k sinh(rj) sinh(r fc )) 2 — 1 
> 0. 

By the same argument as in case 1, we see that d k U, being the graph of the smooth function 
/, is connected and V k , being the region below the positive function / over R> , is also 
connected. □ 

Now back to the proof of proposition 3.2, for part (1), consider the real analytic diffeo- 
morphism u = u(r) : R 3 — > R 3 where Ui = lnr;. The differential 1-form w = X^i^T* 
pulls back (via r = -u _1 (r)) to w — Yli=i @idui as appeared in theorem 3.1. By lemma 
3.3, the extension w = Yli=i @idui is obtained from w by extending each coefficient 6,i by 
constant functions on R 3 — u~ l {U). Thus, by corollary 2.6, the function F(u) = f^w is a 
C 1 -smooth concave function in u G R 3 so that 

(3.7) dF/dm = §i. 

The same argument also works for part (2) since u = u(r) with Ui = lntanh(rj) is a 
real analytic diffeomorphism from R 3 onto R 3 - 



3.3. A proof of theorem 1.4 for Euclidean inversive distance circle packing. 

Suppose otherwise that there exist two inversive circle packing metrics di, c?2 on (S, T) with 
the same inversive distance / G [0, oo)^ so that their discrete curvatures are the same and 
d\ 7^ Ac?2 for any A. Let a G R y be their common discrete curvature. 

We will use the notation that if % G V and x G R y , then Xi = x(i) below. Let be 
the set of all triangles in T. If a triangle s G has vertices k G V, then we denote the 
triangle by s = {i,j, k}. For circle packing metrics of radii r G R^ with a given inversive 
distance /, we use u G H v to denote their logarithm coordinate where Ui = lnr,. Thus, 
there are two points p, q in H v as the logarithmic coordinates of d\ and di so that their 
discrete curvatures are a G R y and p — q ^ A(l, 1, 1, .., 1) for any A. 

We will derive a contradiction by using the locally concave functions F and its concave 
extension F = L w appeared in proposition 3.2 associated to theorem 3.1(1). 

Define a C 1 -smooth function W : H v — > R by 
(3.8) W(u) = - Yl F( Ul ,u,,u k ) + J2( 27r -cn)ui. 

s£T( 2 ),s={i,j,k},i,j,k£V *GV 

The function W is convex since it is a summation of convex functions. Furthermore, 
by the definition of W, (3.7), the definition of discrete curvature (a,), p and q are both 
critical points of W. Since W is convex in H v , p and q are both minimal points of W. 
Furthermore, for all t G [0, 1], tp + (1 — t)q are minimal points of W. In particular, 



W(tp + (l-t)q) = W{p) 
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for all t G [0,1]. Since 



W(tp + (1 - t)q) 



E 



fijk(t) + ^(2tt - ai){tpi + (1 - £)&) 



seT( 2 ),s={i,j,fc},i,i,fcey 



where the function 



(3.9) 



fijk(t) 



F(tpi + (1 - t)q h tpj + (1 - tp k + (1 - t)gfc) 



is convex, it follows that fijk{t) is linear in t G [0, 1] for all triangle s with vertices i,j, k. 
This is due to the simple fact that a summation of a convex function with a strictly 
convex function is strictly convex. By the assumption that p — q ^ c(l, 1, 1) in R^ 
and the surface is connected, there exists a triangle s with vertices i,j,k G V so that 
(Pi,Pj,Pk) ~ (Qi,qj,Qk) ^ (c, c,c) for all c G R. By the given assumption, (pi,Pj,Pk) and 
(<7ij Qji Qk) are m the domain of definition of w in theorem 3.1. Thus for t G [0, 1] close to 
or 1, by theorem 3.1 on the local strictly convexity of —F(ui,v,2,U3) on u\ + ui + = 
and F(u + (c, c, c)) = F(u), fijk{t) is strictly convex in t near 0, 1. This is a contradiction 
to the linearity of fijk(t)- 

3.4. A proof of theorem 1.4 for hyperbolic inversive distance circle packing. 

The proof is essentially the same as in §3.3 and is simpler. For any r G R> > define 
u = u(r) G R^g by U{ = lntanh(rj/2)). For a circle packing with radii r G R>o> let 
u = u(r) and call it the w-coordinate of the circle packing metric. 

We use the same notation as in §3.3. Suppose the result does not hold and let p ^ q G 
R^ be the u-coordinates of the two distinct hyperbolic circle packing metrics having the 
same hyperbolic inversive distance / G R> and the same discrete curvature a = (oj) G R y . 

Define the action functional W on R^ by the same formula (3.8) where F is the concave 
function in proposition 3.2 associated to theorem 3.1(2). Then the same proof goes through 
as in §3.3 since in this case, one of fijk{t) is strictly convex for t near and 1. 

4. 2-dimensional Schlaefli Type Action Functionals and Their Extensions 

The following was proved in [12] . The proof is a straight forward calculation. 

Theorem 4.1. Suppose A is a triangle in the Euclidean plane E 2 ; or the hyperbolic plane 
H 2 ; or the 2-sphere S 2 so that its edge lengths are l\, I2, h and its inner angles are 6\, 82, O3 
where di's angle is opposite to the edge of length ^. Let h G R and let Q be the natural 
domain for length vectors appeared in lemma 2.2. 

(1) For a Euclidean triangle, 
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is a closed 1-form on Q. The integral fti hhh \U)h ^ s ^ oca ^U convex in variable 

i~ h 

u = (u±,U2, U3) where Ui = ln/j for h = and Ui = —\- for h ^ 0. Furthermore, 
I-lh h h) Wh ^ s l° ca tty strictly convex in hypersurface u\ + u 2 + u 3 = 0. 
(2) For a spherical triangle, 

is a closed 1-form on Q. The integral J™ Wh is locally strictly convex in u = 

h 

r/2 



(ui, U2, U3) where Ui = ft /2 sin h 1 (t)dt. 



(3) For a hyperbolic triangle, 

3 

smh h+1 (k) 



^ C /2 sin\t)dt 



is a closed 1-form. 
(4) For a hyperbolic triangle, 

"» = S coth^ft/2) * 

zs a closed 1-form. The integral J™ Wh is locally strictly convex in u = (ui,u 2 ,u 3 ) 
where u, t = ft cot\T h ~ l (t/ 2) dt. 

4.1. Recall that the natural domain Q of the edge length vectors is given by Q = 
{(hi hi h) £ ~&>o\h + lj > hi {hji k} = {1, 2, 3}} for Euclidean and hyperbolic triangles and 
^ = {(hi hi h) £ R>ol^ + (? > hi h + h + h < 2vr, {«, j, k} = {1, 2, 3}}. Let J be the natural 
interval for each individual length /j, i.e., J = R>o for Euclidean and hyperbolic triangles 
and J = (0, 7r) for spherical triangles. In each case of theorem 4.1, there exists a real analytic 
diffeomorphism g : J — >■ g(J) from J onto the open interval g(J) so that Ui = g(k). To 
be more precise, g(t) = hit in the case of /i = of theorem 4.1(1), g(t) = — t — (h 7^ 0) 
in the case of h 7^ in theorem 4.1(1), g(t) = f^ 2 sm~ h ~ 1 (x)dx in the case (2) of theorem 

4.1, g(t) = f* smh~ h ~ 1 (x)dx in the case (3) of theorem 4.1 and g(t) = f* coth~ h ~ 1 (x)dx in 
the case of (4). The real analytic diffeomorphism u(li,l 2 ,h) = (ui,u 2 ,u 3 ) where Ui = g(li) 
sends J 3 onto the open cube g(J) 3 in R 3 . 

By lemma 2.2, each of the angle function 9i(l) : f2 — > R can be extended by constant 
functions to a continuous function : J 3 — > R. Define a continuous 1-form Wh on J 3 by 
replacing 9i in the definition of Wh in theorem 4.1 by 9^. 

Lemma 4.2. The continuous differential 1-form ulh is closed in J 3 . 

Proof. By proposition 2.4 where we take X = J 3 and A = Q, it suffices to show that 
ulh is closed in each connected component U of J 3 — Q. By theorem 4.1 w\a is closed, the 
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restriction of Wh to U is of the form Y^=i °idui where Ui = g(k) and q is a constant. Thus 
Wh\u is closed. □ 

Proposition 4.3. The pull back 1-form on g(J) 3 is a closed 1-form. Further- 

more, if F(u) = f u Wh is locally convex in u(Q) (i.e, in the case (1),(2), (4) of theorem 
4-1), then F(u) = f u {u~ l )*{w h ) is convex in u in g(J) 3 ■ 

Note that by the construction, if u G u(Q) and Wh = Y^i=i a i,h(u)dui (as shown in 
theorem 4.1) then 

(4.i) -asr = «u(«)- 

Furthermore, by definition, the (ph and iph curvatures are sum of two of aj i / l (w)'s. 

PROOF. By corollary 2.6 where we take X = g(J) 3 and A = u(Q), it suffices to show 
that u(Q) is bounded by a real analytic surface in X and F(u) is convex in u(Q) and in 
each component of g(J) 3 — u(Q). 

Since Q in J 3 is bounded by hyperplanes and u(l) = (g(h) , gih) , 9(h)) is a re& l analytic 
diffeomorphism, it follows that u(Q) is bounded by a real analytic surface in g(J) 3 . 

By the assumption F(u) is convex in u(Q). If U is a connected component of g(J) 3 — 
u(fl), then F(u) is linear on U since its partial derivatives are constants on U by the 
construction. Thus by corollary 2.6, the result follows. □ 



5. A Proof of Theorem 1.2 

The argument is essentially the same as that in §3.3. Recall that E is the set of all 
edges in the triangulated surface (S,T). If x e R £ and i e E, we use Xi to denote x(i). If 
s G T^ is a triangle with edges k £ E, we denote it by s = {i,j, k}. 

5.1. A proof of theorem 1.2(3). Suppose otherwise that there exist two distinct 
hyperbolic polyhedral metrics on (S, T) so that their ip h curvatures are the same. Let 
a = (di) G H E be their common iph curvature. 

Recall that a polyhedral metric on (S, T) is given by its edge length map /:£'—>■ R>o- In 
using the variational principle in theorem 4.1(4), the natural variable is given by u : E — > R 
where u(e) = g(l(e)) with g(t) = coth h+1 (s/2)ds. We call it the -u-coordinate of the 
polyhedral metric I and we will use the -u-coordinate to set up the variational principle. 
Therefore there are two distinct points (as -^-coordinates) p ^ q G g(R >0 ) E so that their 
corresponding ip h curvatures are the same a G H E . We will derive a contradiction by using 
the locally strictly convex functions F and its convex extension F introduced in proposition 
4.3 (associated to theorem 4.1(4)). 
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Define a C 1 -smooth function W : g(R > o) E — > R by 

w i u ) = ^2 ^iu h u i-> Uk "> ~ S aiUi - 

seT( 2 ),s={i,j,k},i,j,keE i£E 

The function W is convex since it is a summation of convex functions. Furthermore, 
by the definition of W, (4.1), the definition of ip h and (a,), p and g are both critical points 
of W. Since W is convex, p and q are both minimal points of W. Furthermore, for all 
t G [0, 1], tp + (1 — t)g are minimal points of W. In particular, 

W(tp+(l-t)q) = W(p) 

for all t e [0,1]. Since 

W{tp + (1 - t)g) = Mt) - a ^ + 

i,j,k£E,{i,j,k}£T( 2 ) i&E 

where the function 

(5.1) / ijfe (t) = F(t Pi + (1 - t)%, tp, + (1 - t) qj ,tp k + (1 - 0%) 

is convex, it follows that fijkit) is linear in t e [0, 1]. Since p q, there exists a triangle 
with edges k e E so that (pi,Pj,Pk) ^ (<Zi, 9j, ?fe)- Thus for t e [0, 1] close to or 1, by 
theorem 4.1 on the local strictly convexity, fij k (t) is strictly convex in t near 0, 1. This is 
a contradiction to the linearity of fijkit). 

5.2. A proof of theorem 1.2(2). The proof is exactly the same as above using the 
extended convex function F in proposition 4.3 associated to theorem 4.1(2). 

5.3. A proof of theorem 1.2(1). The proof is the same as that in §5.1 using the 
similarly defined function W. To be more precise, let g{t) = - l — for h ^ and g(t) = Int. 
By the same set up as in §5.1, we conclude that fijkit) given by (5.1) is linear in t. We 
claim this implies that the two Euclidean polyhedral metrics -u _1 (p) and u^iq) differ by a 
scalar multiplication. There are two cases to be discussed depending on h = or h ^ 0. 

CASE 1. h = 0. In this case, p — q ^ c(l, 1, 1) in H E for any constant c. By 
the connectivity of the surface S, there exists a triangle with edges k G E so that 
iPi,Pj,Pk) — (<?i, qj, qk) 7^ (c, c, c) for any constant c. On the other hand, by theorem 4.1(1), 
the action functional F is strictly locally convex in the hyperplane U\ + «2 + «3 = and 
F[u + (c, c, c) = F{u) for a scalar c and u G w(f2). In particular, this implies that the 
function fijkit) is strictly convex in t G [0, 1] for t close to or 1. This contradicts the 
linearity of the function fijkit)- 

CASE 2. h ^ 0. In this case, p ^ cq for any constant c. In particular, there exists 
a triangle with three edges k E E so that iPi,Pj,Pk) 7^ ciq iy qj,qk) for any c G R. By 
theorem 4.1(1), the function fijkit) is strictly convex in t G [0, 1] for t close to or 1. This 
contradicts the linearity of the function fijkit). 

Thus the two polyhedral metrics differ by a scaling. 



RIGIDITY OF POLYHEDRAL SURFACES, III 



17 



References 

[I] Andreev, E. M., Convex polyhedra in Lobachevsky spaces. (Russian) Mat. Sb. (N.S.) 81 (123) 1970 
445-478. 

[2] Bowers, Philip L.; Stephenson, Kenneth, Uniformizing dessins and Belyi (maps via circle packing). 
Mem. Amer. Math. Soc. 170 (2004), no. 805, xii+97 pp. 

[3] Bowers, Philip L.; Hurdal, Monica K. Planar conformal mappings of piecewise flat surfaces. Visualiza- 
tion and mathematics III, 334, Math. Vis., Springer, Berlin, 2003, 

[4] Bobenko, Alexander; Pinkall, Ulrich; Springborn Boris, Discrete conformal maps and ideal hyperbolic 
polyhedra, |http://front. math. ucdavis.edu/10 05. 2698 

[5] Cohn, Henry; Kenyon, Richard; Propp, James, A variational principle for domino tiling. J. Amer. Math. 
Soc. 14 (2001), no. 2, 297-346. 

[6] Chow, Bennett; Luo, Feng, Combinatorial Ricci flows on surfaces. J. Differential Geom. 63 (2003), no. 
1, 97-129. 

[7] Colin de Verdiere, Yves, Un principe variationnel pour les empilements de cercles. Invent. Math. 104 
(1991), no. 3, 655-669. 

[8] Glickenstcin, David, Discrete conformal variations and scalar curvature on piecewise fiat two and three 
dimensional manifolds, to appear in J. of Diff. Geom., |http://front.math.ucdavis.ed u/0906.1560| 
[9] Guo, Ren, Local rigidity of inversive distance circle packing, to appear in Trans AMS, 
|http://front.math.ucdavis.edu/0903.1401| 

[10] Guo, Ren; Luo, Feng, Rigidity of polyhedral surfaces, II, Geom. Topol. 13 (2009), no. 3, 12651312 

[II] Leibon, Gregory, Characterizing the Delaunay decompositions of compact hyperbolic surfaces. Geom. 
Topol. 6 (2002), 361-391. 

[12] Luo, Feng, rigidity of polyhedral surfaces, to appear in J. Differential Geom. 
http://front.math.ucdavis.edu/0612.5714 

[13] Luo, Feng, On Teichmuller spaces of Surfaces with boundary, Duke Journal of Math, 139, no. 3 (2007), 
463-482. 

[14] Luo, Feng, A characterization of spherical polyhedral surfaces. J. Differential Geom. 74 (2006), no. 3, 
407424. 

[15] Rivin, Igor, Euclidean structures on simplicial surfaces and hyperbolic volume. Ann. of Math. (2) 139 
(1994), no. 3, 553-580. 

[16] Stephenson, Kenneth, Introduction to circle packing. The theory of discrete analytic functions. Cam- 
bridge University Press, Cambridge, 2005. 

[17] Thurston, William, Geometry and topology of 3-manifolds, lecture notes, Math Dept., Princeton 
University, 1978, at www.msri.org/publications/books/gt3m/ 

Department of Mathematics, Rutgers University, New Brunswick, New Jersey 08854 
E-mail address: fluo@math.rutgers.edu 



